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The free-energy extrema governing the magnetization-reversal process for a model of an iron 
nanopillar are investigated using the projective dynamics method. Since the time evolution of the 
model is computationally intensive, one could question whether sufficient statistics can be generated, 
with current resources, to evaluate the position of the metastable configuration. Justification for 
the fine-grained discretization of the model that we use here is given, and it is shown that tractable 
results can be obtained for this system on realistic time scales. 
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I. INTRODUCTION 

Magnetic nanopillars can be fabricated^, which offer 
interesting opportunities for direct comparison with nu- 
merical models. In the application discussed here, the pil- 
lars are grown perpendicular to a surface. For an applied 
field parallel to the long axis of the pillar, there exists a 
stable free-energy minimum when the field is parallel to 
the average magnetization of the pillar, and a metastable 
minimum for the antiparallel orientation. Magnetization 
reversal involves a transition from the metastable min- 
imum to the stable minimum across a saddle point in 
the free energy. It is the free-energy difference between 
the saddle point and the metastable minimum that deter- 
mines the switching time of the pillar. While determining 
the magnetization of the metastable minimum is com- 
paratively easy, the magnetization of the saddle point is 
more difficult to determine. Projective Dynamics (PD) is 
a technique that has proven capable of finding the saddle 
point in other magnetization-reversal systems^ 3 -. 



II. MODEL AND NUMERICAL RESULTS 

In order for the proper dynamics to emerge, the mag- 
netization of the simulated pillar must be discretized on a 
sufficiently fine scale. Pillar models that have previously 
been studied by projective dynamics consisted only of a 
one-dimensional chain of spins 3 -. This model lacks richer 
dynamics, such as magnetization curling modes, that are 
seen in experiments and in more realistic models. In ad- 
dition, the dependence of the switching field on the angle 
of misalignment between the pillar axis and the applied 
field observed in the one-dimensional model 4 does not 
correspond to that seen in the experimental pillars***^. 

To achieve more realistic dynamics, a model pillar with 
physical dimensions 10 nm x 10 nm x 150 nm was dis- 
cretized onto a 6 x 6 x 90 cubic lattice. The model 
parameters were chosen to be consistent with bulk iron: 



exchange length of 2.6 nm and magnetization density of 
1700 emu/cm 3 . The temperature was fixed at 20 K. 

The dynamics of the system are governed by the 
Landau-Lifshitz-Gilbert (LLG) equation^, 
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which describes the time evolution of the magnetization 
in the presence of a local field. Here 70 is the electron 
gyromagnetic ratio with a value of 1.67 x 10 7 Hz/Oe, and 
a is a phenomenological damping parameter, chosen as 
0.1 to give underdamped behavior. M[fi) is the mag- 
netization density at position fl, with a constant mag- 
nitude AT,, and H(fl) is the total local field at fl. The 
latter contains contributions from the dipole-dipole in- 
teractions, exchange interactions, and the applied field. 

A reversal simulation begins with the magnetization 
oriented along the long axis of the pillar in the +z di- 
rection. It is allowed to relax in the presence of a field 
+Boz before the field is quickly switched to —Bgz, where 
Bq = 1125 Oe has been chosen smaller than the coercive 
field. 

Projective dynamics involves analyzing the probabil- 
ity of growing toward the stable equilibrium or shrink- 
ing away from it, as a function of a single variable that 
measures the progress of the reversal process. To collect 
growing and shrinking statistics here, the z-component of 
the total magnetization of the pillar, M z , was recorded 
at each integration step. In addition, the change in this 
value from the previous integration step, AM Zl was also 
recorded. The magnetization axis was discretized into 
bins into which the collected data were sorted. The bin 
size was determined such that the discretized AM 2 would 
only involve changes between adjacent bins. The growing 
and shrinking probabilities were calculated by counting 
the number of times that AM Z would cause the magne- 
tization to change bins. Define G as the number of times 
the magnetization hopped to the next bin in the direc- 
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FIG. 1: Growing and shrinking probabilities versus the mag- 
netization along the z-axis for a 6x6x90 pillar. Grow/shrink 
data were collected into 6000 bins, which span the entire mag- 
netization axis, [-1,1]. The simulation was run at 20 K, 1125 
Oe, and the applied field was directed along the long axis 
of the pillar, antiparallel to the initial magnetization. The 
continuous curves represent nine-point running averages 

tion toward equilibrium, and S as the number of hops in 
the direction back toward the metastable state. If N is 
the total number of visits to the starting bin, then the 
growing and shrinking probabilities for each bin, Pq and 
Ps, can be written respectively as 



The crossings of the growing and shrinking prob- 
abilities give information about the locations of the 
metastable free-energy minimum and the saddle point. 
Where the two probabilities are equal, the system has no 
preferred direction of magnetization evolution, indicating 
that the free-energy landscape is flat in this region. This 
is the location of an extremum. The right-hand intersec- 
tion in Fig. H corresponds to the metastable free-energy 
minimum, where the system spends most of its time. At 
the saddle point, the growing and shrinking probabili- 
ties are also equal, as indicated by the left-hand crossing. 
The region between these two crossings shows a higher 
probability for shrinkage than growth, corresponding a 
preference for moving toward the metastable minimum. 

Previous simulations performed on a one-dimensional 
chain of seventeen spins exhibited Stoner-Wohlfarth be- 
havior for the angular dependence of the switching field 4 . 
This is not representative of experimental data£*§*£. How- 
ever, a full three-dimensional model qualitatively pro- 
duced the proper angular dependence 4 . This is not sur- 
prising, given that the full three-dimensional model al- 



lows for magnetization curling modes which are simply 
not allowed in the one-dimensional model. PD has pre- 
viously been applied to the one-dimensional spin chain 3 , 
but several thousand switches had to be performed in or- 
der to accumulate sufficient statistics to reveal the free- 
energy extrema 3 ". With typical switches taking around 
1500 CPU-hours for the full three-dimensional model, it 
was not clear that sufficient statistics could be generated. 

Figure ^ shows the results collected from only 20 
switches for the full model. The estimated Pq is shown 
as circles, and Ps as crosses. Scatter in these data is 
primarily due to the counting statistics and the small 
number of samples in each bin. Nine-point running aver- 
ages are used to reduce the scatter and improve estimates 
of the crossings of Pg and Ps. The running averages are 
represented by curves and are sufficient to determine the 
locations of the extrema with an uncertainty of about 
1 x 10~ 3 in the magnetization. For the one-dimensional 
model the estimated curves representing Pq and Ps are 
almost parallel near the saddle point, making it difficult 
to locate crossings of the growing and shrinking proba- 
bilities (see Fig. 1 of Reff 3 "). Since they cross at a sig- 
nificant angle for the model studied here (see Fig. QJ, it 
is possible to find an estimate for the magnetization of 
the free-energy saddle point from far fewer switches than 
were needed for the one-dimensional model. 



III. SUMMARY AND CONCLUSIONS 

The projective dynamics method was used to produce 
growing and shrinking probabilities during the magne- 
tization reversal of a realistic model nanopillar. These 
probabilities can be used to locate the magnetization of 
the free-energy extrema. In particular, the magnetiza- 
tion of the saddle point is revealed. Compared to the 
one-dimensional model, far fewer switches are needed to 
find the location of the extrema in the three-dimensional 
model. In the future, a study of the dependence of the 
saddle point on temperature and misalignment between 
the pillar and the applied field are planned. 
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